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Functional magnetic resonance imaging (fMRI) techniques allow 
definition of cortical nodes that are presumed to be components of 
large-scale distributed brain networks involved in cognitive 
processes. However, very few investigations examine whether 
such functionally defined areas are in fact structurally connected. 
Here, we used combined fMRI and diffusion MRI-based tractog- 
raphy to define the cortical network involved in saccadic eye 
movement control in humans. The results of this multimodal 
imaging approach demonstrate white matter pathways connecting 
the frontal eye fields and supplementary eye fields, consistent with 
the known connectivity of these regions in macaque monkeys. 
Importantly, however, these connections appeared to be more 
prominent in the right hemisphere of humans. In addition, there was 
evidence of a dorsal frontoparietal pathway connecting the frontal 
eye field and the inferior parietal lobe, also right hemisphere 
dominant, consistent with specialization of the right hemisphere for 
directed attention in humans. These findings demonstrate the utility 
and potential of using multimodal imaging techniques to define 
large-scale distributed brain networks, including those that 
demonstrate known hemispheric asymmetries in humans. 
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Introduction 

It is now widely accepted that many cognitive functions are 
underpinned by large-scale brain networks connecting distrib- 
uted cortical regions, some of which are highly lateralized to 
the left or right hemisphere in humans (Mesulam 1981; 
Goldman-Rakic 1988; Catani and ffytche 2005). Establishing 
how the cortical nodes within such networks are intercon- 
nected is likely to be crucial for understanding their functional 
significance. However, until the recent development of 
advanced imaging methods, we have gained little direct 
knowledge on network connectivity within the human brain. 

Instead, we have largely relied on information from 
anatomical tracer studies in nonhuman primates to infer 
patterns of connectivity in the human brain. Indeed, white 
matter tracts that have been described in detail in macaque 
monkeys (Schmahmann and Pandya 2006) are often cited as if 
they are known to exist in humans, but for the main, have yet 
to be confirmed. Until recently, the only way to identify white 
matter anatomy in humans was through postmortem dissection 
and histological analysis of white matter fibers. However, with 
recent advances in diffusion magnetic resonance imaging 
(MRI) tractography, it is now possible to perform "virtual 
dissections" of the living human brain to visualize large-scale 
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connectivity maps (Mori et al. 1999; Basser et al. 2000; Catani 
et al. 2002; Sporns et al. 2005; Hagmann et al. 2008; Jones 
2008). Yet, very few investigations have examined how cortical 
regions defined on the basis of their functional activation are 
structurally connected. 

Here, we combined tractography with functional MRI (fMRI) 
to investigate structural connections between nodes in the 
cortical network for gaze control in humans. This system is of 
fundamental importance as oculomotor responses can be used 
to probe many aspects of cognitive control (Kennard et al. 
2005; Sweeney et al. 2007) including spatial working memory 
(Husain et al. 2001; Brignani et al. 2010), visual search and 
attention (Binello et al. 1995; Mannan et al. 2005) as well as top- 
down control related to response initiation and suppression 
(Sumner et al. 2006; Anderson et al. 2008). It is now well 
established in nonhuman primates that a distributed network, 
including the frontal eye field (FEF) in the dorsolateral frontal 
cortex (Bruce et al. 1985; Huerta et al. 1987), supplementary 
eye field (SEF) located in dorsomedial frontal cortex (Russo and 
Bruce 2000), and parietal eye field (PEF) within the lateral 
intraparietal (LIP) area (Thier and Andersen 1998; Goldberg 
et al. 2002), is crucial for saccadic eye movement control (for 
a review, see Johnston and Everling 2008). Moreover, anatom- 
ical tracer studies have demonstrated that these areas are 
heavily interconnected, both within and across hemispheres 
(Stanton et al. 1993, 1995; Schall et al. 1995). 

The human homologues of these regions have been 
extensively documented using high-resolution neuroimaging 
techniques (Paus 1996; Grosbras et al. 1999; Tehovnik et al. 
2000; Lobel et al. 2001; McDowell et al. 2008; Amiez and 
Petrides 2009), but despite substantial literature on their 
location and function, the underlying pattern of connections 
remains to be established. Investigations of oculomotor control 
in clinical populations have also highlighted a number of 
deficits that might, in part, be ascribed to pathological changes 
in white matter pathways (Broerse et al. 2001; Gooding and 
Basso 2008) and observations in patients with brain lesions 
have suggested right hemisphere specialization for directed 
attention, a process that is intimately linked to gaze control 
(Posner et al. 1984; Rafal 1994; Mesulam 1999). Hence, 
mapping the structural connections that exist within the 
healthy human brain will help our understanding of the 
changes in circuitry that might explain such deficits. 

Typically, previous attempts to infer structure -function 
relationships in the human brain have used tractography to 
identify pathways connecting anatomically defined regions, 
followed by post hoc interpretation of their functional 
significance based on previously published functional data. 
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Alternatively, fMRI has been used to identify cortical areas 
associated with a specific cognitive function and then 
assumptions about the underlying structural connections have 
been made. However, by combining tractography and fMRI 
within the same group of subjects, it should be possible to 
identify "maps of connections" implicated in a specific cogni- 
tive function (Rykhlevskaia et al. 2008). 

To date, there are very few multimodal imaging studies that 
have combined these 2 techniques within the same subjects 
(e.g., Kim and Kim 2005; Saur et al. 2008; Staempfli et al. 2008; 
Lanyon et al. 2009), and these investigations have typically used 
the single approach of constraining fiber reconstruction 
algorithms to dissect pathways connecting functionally defined 
regions. Here, to advance on these methods, we combined 
tractography and fMRI in 10 healthy volunteers, using a 2-stage 
analysis. For each individual, we initially ran an anatomically 
guided tract reconstruction (independent of the fMRI data), 
which ensured no loss of information by taking into account 
the variability that exists in individual anatomy. This was 
followed by a second stage in which individual fMRI data were 
used to constrain the tractography results. In addition, we 
created a group -averaged diffusion MRI data set to obtain 
a group-averaged tractography result to further explore the 
anatomy that might generalize to the population. 

Materials and Methods 

Participants 

Ten participants (4 males), all right handed (mean age 29 ± 3 years), 
with normal visual acuity, and no history of neurological or 
psychological disorder, gave informed written consent to take part in 
this study. This study was approved by Hammersmith Ethics Commit- 
tee. 

Each participant underwent 4 consecutive scans in a single session, 
using a GE 1.5-T Signa NV/i LX MRI system (General Electric, 
Milwaukee, WI), with standard quadrature birdcage head coil: 1) 
diffusion-weighted echo planar imaging (EPI), 2) functional EPI — eye- 
field localizer A, 3) functional EPI — eye-field localizer B, and 4) T r 
weighted structural scan. 

Diffusion Tensor-MRI Acquisition 

A multislice diffusion-weighted EPI sequence, fully optimized for 
diffusion tensor imaging of white matter (Jones, Williams, et al. 
2002), was carried out on all subjects. Full brain coverage was achieved 
with 60 near-axial slices, with isotropic resolution (voxel size 2.5 x 2.5 
x 2.5 mm, field of view 240 x 240 mm, acquisition matrix 96 x 96, 
interleaved slice order). The duration of the diffusion-encoding 
gradients was 17.3 ms, giving a maximum diffusion weighting of 1300 
s mm" 2 . Diffusion gradients were applied in 64 isotropically distributed 
orientations. The acquisition was gated to the cardiac cycle using 
a peripheral gating device placed on the subjects forefinger. Depending 
on the individual heart rate, an effective TR was selected that allowed at 
least 180 ms per image, that is, enough time to collect, reconstruct, and 
save each slice. Scan time was -20 min, depending on heart rate. 

Analysis of Diffusion Tensor-MRI data 

After correction for the image distortions introduced by the application 
of the diffusion-encoding gradients (Leemans and Jones 2009), the 
diffusion tensor was determined in each voxel following the method of 
Basser et al. (1994). From the diffusion tensor data, information about 
the average diffusivity, diffusion anisotropy, and the principle direction 
of diffusion can be calculated (Basser and Pierpaoli 1996; Basser et al. 
2000). The Diffusion Tensor-MRI (DT-MRI) data were coregistered to 
each individual's structural scan and registered in Montreal Neurolog- 
ical Institute (MNI) stereotactic space for consistency with the fMRI 



data. The fractional anisotropy (FA) was computed in each voxel, 
highlighting the boundary between white and gray matter. This "FA 
map" was used to define the regions of interest (ROIs) used in the 
tractography process. 

To ensure that no hemispheric bias could have been imposed by the 
observer whilst defining the ROIs, half of the diffusion MRI data sets 
were flipped along the midline (as well as their corresponding 
anatomical and fMRI data) before the ROIs were defined. Information 
about which data sets were flipped was not revealed to the observer 
until after the analysis was complete, at which point all data sets were 
returned to their original orientation for statistical analysis and 
presentation of the results. 



Group Averaged DT-MRI Data Set 

A group-averaged DT-MRI data set was generated from all 10 
participants. The (unflipped) tensor data for each participant was 
spatially normalized and coregistered to the standard T 2 -weighted MNI 
EPI template, included as part of SPM99 (Wellcome Trust Centre for 
Neuroimaging, UCL, www.fil.ion.ucl.ac.uk), using the approach de- 
scribed by Alexander et al. 1999, which employs the Automated 
Image Registration (AIR) package for coregistration (Woods, Grafton, 
Holmes, et al. 1998; Woods, Grafton, Watson, et al. 1998). Full details of 
the spatial normalization and coregistration procedure and generation 
of the average DT-MRI volume are described in Jones, Griffin, et al. 
(2002). 



Eye-Field localizer and JMRI Acquisition 

All subjects underwent 2 functional EPI scans to localize the cortical 
eye fields. Subjects were required to perform blocks of large self-paced 
horizontal saccades in the dark for 30 s, followed by 30 s of rest, in 
which they maintained a steady central eye position. A practice session 
prior to scanning was completed by all subjects to ensure they could 
perform the task correctly. The location of the cortical eye fields differs 
according to task requirements, stimulus type, and the oculomotor 
response required (McDowell et al. 2008). For these reasons, we used 
an oculomotor paradigm that minimized target-related effects and 
hence minimized activity unrelated to generating and carrying out 
a saccade. Aural instructions were given via MRI compatible head- 
phones instructing participants when to "start" and "stop" making eye 
movements. 

T 2 -weighted images were acquired using Gradient Echo EPI (128 x 
128 matrix, FOV 240 x 240 mm, time echo [TE] 40 ms and 90° flip 
angle). Full brain coverage was achieved with 28 near-axial slices (4.5- 
mm slice thickness with a 0.5 -mm gap, interleaved slice order, in-plane 
resolution 1.9 x 1.9 mm). Participants performed 2 identical scan runs, 
each lasting 8 min (96 volumes per scan, TR 5 s). Each scan started with 
a 30 s "rest" block (6 volumes), followed by a 30 s "eye movement" 
block (6 volumes), repeated 8 times. 

A T x -weighted anatomical scan with 1.5-mm slice thickness and 0.86 
x 0.86 in-plane resolution (after zero filling) was also acquired for every 
participant using an axial 3D SPGR sequence (256 x 256 matrix, FOV 
220 x 176, TR 159 ms, TE 5.2 ms, 20° flip angle). 

Analysis of JMRI Data 

All imaging data were analyzed using SPM2 (www.fil.ion.ucl.ac.uk). For 
the individual analyses, half of the fMRI data sets were flipped along the 
midline prior to preprocessing — corresponding to those participants 
who had their DT-MRI data sets flipped. The corresponding structural 
scan was also flipped. This was to ensure that no hemisphere bias was 
imposed by the observer when defining the ROIs. 

For each subject, all images were realigned to the first image to 
compensate for head movement, coregistered to the participant's T x 
structural image, and spatially smoothed with a 7-mm isotropic 
Gaussian smoothing kernel. For the group analysis (no data sets were 
flipped), all images for each subject were realigned to the first image, 
spatially normalized to the SPM EPI template, and spatially smoothed 
with a 7-mm isotropic Gaussian smoothing kernel. 

For each individual, linear regressors representing the time series for 
blocks of eye movements and blocks of rest were convolved with 
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a synthetic hemodynamic response function and its temporal de- 
rivative. The general linear model, as employed by SPM2, was used to 
generate parameter estimates of activity for each regressor at every 
voxel. A linear contrast between regressors identified those voxels with 
greater activity during blocks of eye movements compared with rest, 
thus revealing the locus of the cortical eye fields in each subject. A 
threshold of P < 0.001 uncorrected was used to determine significance. 

For the group data, a random effects analysis was performed (Friston 
et al. 1999). A single mean image representing "eye movements versus 
rest" for each subject was used as the basis for interparticipant 
comparisons and used to generate a statistical parametric map of the t 
statistic at every voxel (Fig. 1). A threshold of P < 0.001 uncorrected 
was used to determine significance. 



Anatomically Defined ROIs 

For stage 1 of the analysis, all ROIs were defined on axial slices of the FA 
map, which clearly highlights the boundary between white and gray 
matter (Fig. 2). Initially, for all subjects, 2 large ROIs were defined in the 
superior frontal lobe (SEF ROI) and middle frontal lobe (FEF ROI) of 
each hemisphere, using anatomical landmarks for guidance. ROIs were 
defined to encompass a large region of white matter lying directly 
beneath regions of cortex known to include the human SEF and FEF. 
For the FEF ROI, a region of white matter underlying the middle and 
superior frontal gyrus, lateral to the internal capsule, and anterior to the 
central sulcus, was defined (approximate Brodmann areas [BAs] 44, 8, 6, 
and 4). The approximate MNI coordinates for this region ranged in the 
X direction from 25 to 60, in the F direction from -10 to 25, and from 
20 to 60 in the Z direction. These coordinates varied depending on 
individual anatomy. For the SEF region, a region of superior frontal 
cortex medial to the internal capsule and anterior to the central sulcus, 
underlying the precentral gyrus (BAs 4 and 6), was defined 
(approximate X range: 4-20, Frange: -15 to 10, Z range: 45-70). 
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Figure 1. Functionally defined cortical eye fields. The cortical eye fields were 
localized in each participant using a standard saccadic eye movement paradigm (see 
Materials and Methods). Group normalized data are presented here, illustrating 
regions of the brain where activity was greater during blocks of saccadic eye 
movements in the dark compared with a central fixation baseline. Random effects 
analysis at P = 0.001 was used to determine coordinates of activation maxima, Z 
scores and approximate BAs as given above. A threshold of P = 0.005 has been 
used for illustration purposes only. 



A further ROI was then denned within the parietal lobe of each 
hemisphere, encompassing a large area of white matter underlying all 
regions of the intraparietal sulcus (IPS) known to be involved in 
saccadic eye movements — that includes IPS 1, 2, and 3 (Sereno et al. 
2001; Schluppeck et al. 2005, 2006; Levy et al. 2007; Anderson et al. 
2008). More specifically, the PEF ROI started at the level of the corpus 
callosum inferiorly and included all white matter posterior to the 
occipitotemporal border (where regions BA 19 and 39 meet, MNI 
-60), extending from the cuneaus on the medial side to the angular 
gyrus laterally and from the middle occipital gyrus inferiorly to the 
inferior aspect of the IPS as well as extending anteriorly along the IPS 
(approximate X range: 0-45, Frange: -50 to -90, Z range: 30-60). 

For the group analysis, similar anatomical ROIs were defined on the 
group-averaged FA map and used as seed points for running 
tractography on the group-averaged DT-MRI volume. 



Functionally Redefined ROIs 

For stage 2 of our analysis, we redefined our regions using each 
individual's fMRI data. Using visual markers, the original large 
anatomically defined ROIs were refined to only include regions of 
white matter that lay directly beneath areas of cortical activity. These 
ROIs were smaller in si2e and more closely respected the anatomical 
locality of the individual's eye fields. The same tractography algorithms 
were then rerun for these new ROIs and the reconstructed streamtube 
data visualized in 3D along with the fMRI data. From this interactive 3D 
visualization, we were able to confirm whether the dissected pathways 
accurately reached areas of cortical activity. 

Similar to the individual analyses, the large anatomically defined ROIs 
used for the group analysis were refined using the group normalized 
fMRI results and tractography rerun on the group-averaged DT-MRI 
volume. 



Tractography Algorithm 

White matter fiber trajectories were reconstructed and visualized in 3D 
using in-house software based on the procedure originally described by 
Basser et al. (2000). Whole brain tractography was performed, by 
seeding tracts from every vertex of a 2 x 2 x 2 mm grid superimposed 
on the brain. For each of these seed points, the principal eigenvector 
was calculated from the DT data (i.e., the direction of greatest 
diffusivity), which was treated as being tangential to the trajectory of 
the tract. Each fiber trajectory was reconstructed by propagating 
a streamline bidirectionally from an initial seed point in the direction of 
the principal eigenvector for 0.5 mm and then recomputing the 
principal eigenvector at the next location. The algorithm then moved 
a further 0.5 mm along this new direction. The tracking process 
continued tracing a pathway in this manner until the FA of the tensor 
fell below a fixed arbitrary threshold (set to 0.2). This threshold 
consistently differentiated between white and gray matter when 
windowing the FA image. 

Only those trajectories that passed through 2 predefined ROIs, 
defined on the FA images, were retained for analysis. Using this 2 -ROI 
approach, tractography was used to reconstruct all fibers passing 
through any 2 of our spatially distinct ROIs. In particular, we were 
interested in dissecting tracts which connect 1) SEF and FEF, 2) FEF 
and PEF, and 3) SEF and PEF. The dissected tracts were visualized in 
3D as "streamtubes," using the mathematical software package 
MATLAB (The Mathworks, Natick, MA) (Zhang et al. 2003). This 3D 
visualization could also be displayed simultaneously with the FA 
volumes, sliced into axial, coronal, and sagittal planes, with or without 
superimposed fMRI activity and viewed interactively from any angle 
to facilitate identifying the neuroanatomical location of the tract 
reconstructions. 

For each dissected pathway, we computed the number of recon- 
structed streamlines (RS) and the mean FA. 

Results 

First, we obtained diffusion-weighted EPI brain scans on 10 
healthy volunteers (all right handed), using a sequence fully 



Cerebral Cortex April 2012, V 22 N 4 767 



optimized for diffusion MRI tractography Qones et al. 1999). In 
addition, we functionally identified the cortical eye fields in 
each of our 10 subjects using functional MRI and a standard eye 
movement paradigm (see Materials and Methods) (Fig. 2). One 
subject had very little cortical activity above threshold in any 
brain area and was therefore discarded from the individual 
analysis. For the remaining 9 subjects, we carried out a 2-stage 
process for dissecting the white matter tracts connecting the 
cortical eye fields. 



Stage 1: Diffusion MRI- Based Tractography 
Reconstruction of the Connections between Anatomically 
Defined Frontal Oculomotor Regions 

In this analysis, we reconstructed the white matter pathways 
of interest independently from the fMRI results. By using only 
anatomical landmarks, we were able to visualize frontal 
connections both within and across hemispheres. This 
analysis allowed us to explore the individual anatomy of the 
white matter pathways and therefore assess interindividual 
variability in patterns of connections. We first defined 2 large 
ROIs within the dorsolateral and dorsomedial prefrontal 
cortex of each hemisphere for all subjects. These regions 
were defined on the individual FA maps, using anatomical 
landmarks for guidance and were drawn to encompass white 
matter underlying regions of cortex known to include the 
FEF and the SEF. For full details of the location and extent of 
the ROIs defined, see Materials and Methods. Whole brain 
tractography was then performed (seeding tracts from every 
vertex of a 2 x 2 x 2 mm grid superimposed on the brain), 
but only those trajectories that passed through both ROIs 



were retained for analysis (Conturo et al. 1999; Catani et al. 
2002). In this way, we dissected the white matter pathways 
between the FEF and SEF within each hemisphere, as well as 
the pathway between the right and left SEF across 
hemispheres. 

Comparing the number of voxels within each ROI for the 
right and left hemisphere confirmed there was no significant 
difference in region size between hemispheres for the FEF 
(fe= 0.527, P = 0.613) or SEF (fe = -0.265, P = 0.798). Using 
these ROIs, we were able to dissect white matter pathways 
directly connecting the dorsomedial ROI (SEF) and the 
dorsolateral ROI (FEF) within each hemisphere for all 9 
subjects, as well as a pathway connecting the right and left 
SEF via the anterior body of the corpus callosum. The 
dissected pathways were viewed in 3D, simultaneously with 
each individual's fMRI data and superimposed on intersecting 
axial and coronal slices of the FA map — allowing the origin, 
course, and termination of each pathway to be determined 
(Fig. 3A). For all 9 subjects, the path between the SEF and FEF 
in both hemispheres accurately terminated in regions of 
cortical activity, as did the pathway connecting the right and 
left SEF across hemispheres. Interestingly, the number of 
RS in this pathway was significantly higher in the right 
hemisphere than the left (t(8)= 2.770, p = 0.024) (Fig. 4). 
Note — this is despite there being no hemispheric bias in the 
size of the ROIs (as reported above). 

In addition to the reconstructed pathways that accurately 
terminated within regions of cortical activity, there were also 
a number of pathways that projected to adjacent areas. This is 
not unexpected given the large size of our initial ROIs. To 
ensure that the hemispheric bias remained reliable and robust 






ROI 



Figure 2. Anatomical definition of the ROIs. For stage 1 of the analysis, all ROIs were defined on the individual FA maps, using anatomical landmarks for guidance and were 
drawn to encompass white matter underlying regions of cortex known to include the FEF, SEF, and PEF. Selected axial slices of the FA map for a representative subject are 
illustrated above, with voxels to be included in the right FEF, SEF, and PEF ROIs outlined in orange. Note that the frontal and parietal ROIs do not cross the central sulcus (red line) 
and the FEF and SEF ROIs do not cross the internal capsule. For clarity, not all slices are illustrated in this figure — for full details of theX, Y, and Z ranges included in each ROI, see 
Materials and Methods. The Z coordinates given are approximate, as the individual FA maps were not normalized to the MNI template for the individual subject analyses. 
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Figure 3. Combined DT tractography and fMRI for a single subject. For a representative subject, reconstructed tracts are illustrated in 3D superimposed on coronal, axial, and 
sagittal slices of the individual's anatomical image, together with their fMRI data (illustrated at P = 0.005). (4) Connections within the frontal lobes — contralateral SEF to SEF 
tract illustrated in red, ipsilateral SEF to FEF tracts in green. [B) Lateral frontoparietal pathway — dorsolateral tract connecting the right FEF and inferior parietal cortex in yellow. (C) 
Medial parietal connections — tract connecting ipsilateral SEF and PEF in magenta and contralateral PEF to PEF connection in green. This figure was generated using "ExploreDTI" 
(Leemans et al. 2009). 
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Figure 4. Histogram showing the number of reconstructed streamtubes in each tract. Group mean number of streamtubes in each reconstructed tract for {A) anatomically 
defined seed regions, (B) functionally redefined seed regions (see Materials and Methods). Error bars indicate standard error of the mean. 



for only those pathways that accurately terminated within the 
functionally defined FEF and SEF, we performed a second stage 
of analysis in which our original ROIs were redefined using 
each individual's fMRI activity for guidance. We then reran the 
tractography using these new ROIs. 

Stage 2: Tractography Constrained Using Frontal JMRI 
Activity to Redefine Seed Points 

For the second stage of analysis, we refined the anatomical 
ROIs to only include voxels of white matter lying directly 
beneath the functionally defined FEF and SEF for each 
individual (a threshold of P < 0.001 uncorrected was applied 
to each individuals fMRI data to determine significance). 
There was no significant difference in ROI volume between 
the new right and left FEF ROI (7 8 = 0.478, P= 0.647) or SEF 
ROI (fe= 0.049, P= 0.962). 

Using these new functionally defined ROIs, we reran our 
tractography analysis. The original reconstructed pathways 
were all found to be reliable, and although there were now 
fewer RS within each pathway, the right hemisphere lateral- 
ization remained significant for the SEF to FEF connection 
(Fig. 4B) = 3.407, P= 0.011). Furthermore, for each subject, 
we obtained the mean value of the FA sampled along each 
pathway. This is a scalar measure that reflects the degree to 
which the diffusivity depends on the orientation in which it is 
measured and therefore considered an index reflecting 
microstructural organization (e.g., cohesiveness, ordering, 
etc.) and biological properties (e.g., degree of myelination, 
membrane permeability, etc.) of fibers (Beaulieu 2002). 
Although there were no significant differences between the 2 
hemispheres at a threshold of P = 0.05, there was a trend for 
the mean FA to be larger on the right than the left for the SEF to 
FEF tract (fe = 1.967, P= 0.085). 

Anatomical and JMRI- Constrained DT Tractography of 
the Frontoparietal Connections 

Posterior projections from the FEF to the LIP area in the IPS 
have also been identified in nonhuman primates (Huerta et al. 
1987; Blatt et al. 1990; Schall et al. 1995; Stanton et al. 1995). 



Hence, having successfully reconstructed white matter path- 
ways connecting the 2 eye fields of the frontal cortex in our 9 
subjects, we then sought to determine whether posterior 
projecting pathways could also be established. 

The human homologue of monkey LIP is thought to lie within 
the IPS (Sereno et al. 2001; Schluppeck et al. 2005, 2006; Levy 
et al. 2007), where topographically organized areas associated 
with saccadic eye movements have been identified. Therefore, 
using the same methods as before, we initially defined a large 
parietal ROI in the white matter underlying the proposed location 
of the human PEFs in both hemispheres. Statistical analysis 
confirmed there was no significant difference in the size of the 
right and left parietal ROIs (fe = -0.231, P= 0.823). Using the same 
2 -ROI approach, we attempted to dissect white matter tracts 
connecting the FEF to PEF and the SEF to PEF within each 
hemisphere, as well as a tract connecting the right and left PEF 
across hemispheres. 

We were able to reconstruct a dorsal frontoparietal tract 
connecting the FEF ROI to a posterior parietal region in both 
hemispheres, for all 9 subjects (Fig. 3E). This pathway 
consistently had a greater number of RS in the right hemi- 
sphere compared with the left (^ = 3.566, P= 0.007) (Fig. 4 A). 
However, examining the reconstructed pathway in 3D along 
with the fMRI data revealed that this pathway did not project to 
a region of parietal cortex consistently activated by our 
oculomotor task (when a standard threshold for significance 
was applied: P < 0.001). Instead this pathway terminated in 
a more lateral region of cortex, known to be involved in visually 
guided oculomotor tasks (Perry and Zeki 2000; Mort, Perry, 
et al. 2003; Petit et al. 2009). The parietal ROIs used as seed 
points for the tractography were defined anatomically to 
include white matter underlying all regions of the IPS known 
to be involved in saccadic eye movements — that includes IPS 1, 
2, and 3 (Sereno et al. 2001; Schluppeck et al. 2005, 2006; Levy 
et al. 2007; Anderson et al. 2008). Hence, this ROI extended 
from posterior to anterior ends of the IPS, thereby incorporat- 
ing white matter associated with this dorsolateral frontopar- 
ietal tract (see Discussion). 
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A dorsomedial pathway connecting the SEF and PEF ROIs, via 
the cingulum, was also found for all subjects in the left 
hemisphere, but for only 6 subjects in the right hemisphere, 
although there was no consistent hemisphere bias in the 
number of streamtubes (fe = 0.645, P= 0.540) (Figs. 3Cand 4B). 
Finally, as expected, there was a U-shaped tract connecting the 
contralateral PEF via the posterior body of the corpus callosum 
(Fig. 36). All of these tracts accurately terminated in regions of 
cortical activity. 

When the tractography was rerun using the more conser- 
vative ROIs denned by fMRI data, only -50% of fibers within the 
large pathway connecting the FEF to lateral parietal cortex 
survived, and the rightward hemisphere bias was abolished (t s = 
0.635, P = 0.546) (Fig. 4B). In contrast, the more variable 
pathway connecting the SEF and PEF did survive, but due to 
this pathway not being reconstructed in all subjects, we cannot 
make reliable statistical inferences. Furthermore, there was no 
significant hemisphere bias in the mean FA for either the FEF to 
PEF pathway (fe = -1.208, P= 0.272) or the SEF to PEF pathway 
(fe = -1.207, P= 0.294). 

Group Average Analysis of the Diffusion and JMRI data 

Lastly, we created a group average diffusion MRI data set to 
obtain a group-averaged tractography result (Jones, Griffin, et al. 
2002). Group-averaged tractography has the advantage of 
providing a smooth and continuous trajectory representing the 
"central" location of the white matter pathways within the group, 
from which we can derive assumptions about the connections of 
the cortical eye fields within the general population. A further 
advantage of the group analysis is that the data are normalized 
and coregistered to a common space (MNI space, which 
approximates to the space of Talaraich and Tournoux), which 
allows for easy comparison with previously published coordi- 
nates for the location of the human cortical eye fields, as well as 
allowing our reconstructed pathways to be assigned to currently 
established white matter pathways. 

We reran the tractography using the same 2-ROI approach as 
before, to obtain a group-averaged tractography result for 
connections between the SEF, FEF, and PEF. As before, the 
reconstructed tracts were visualized in 3D together with the 
group -averaged fMRI results (Fig. 5). Consistent with the 
individual data, within the frontal lobes there was a direct 
pathway connecting the ipsilateral FEF and SEF which was right 
hemisphere lateralized (number of streamtubes: right-.left = 
2.60:1). There was also a direct tract connecting the 
contralateral SEF via the corpus callosum (Fig. 5v4). Posterior 
projections from the FEF were similar to those obtained for the 
large anatomically defined ROIs in the individual data sets, with 
a dorsal pathway connecting the FEF and a lateral area of 
parietal cortex. Similar to the individual data, this tract was also 
right hemisphere lateralized (number of streamtubes: right:left 
= 2.55:1), and once again did not appear to terminate in 
a region of cortex activated by our eye movement localizer 
(coordinates of tract end point: R and L [±39/-58/46], 
homologous BA 39/40) (Fig. 5B) but, is known to be involved 
in visually guided saccades. Hence, this right lateralized tract, 
which appears to generalize to the population, may also play 
a role in oculomotor control. For the group data, we were 
not able to reconstruct the medial tract connecting the SEF 
and PEF that had been evident in a subset of our individual data 
sets. 



Discussion 

By combining fMRI and diffusion MRI-based tractography, we 
have reconstructed the candidate white matter pathways 
linking functionally specified oculomotor areas in the human 
brain. Our data confirms some degree of homology with white 
matter connections known to exist in macaque monkeys but 
also highlights important characteristics unique to humans, 
particularly with respect to hemispheric differences. 

For all individuals, we reconstructed pathways between the 
right and left SEF that passed through the anterior body of the 
corpus callosum (Figs 3 A and 5 A) and the right and left PEF via 
the posterior body of the corpus callosum (Figs 3C and 5B). 
There was also a direct pathway connecting the SEF to the 
ipsilateral FEF in both hemispheres, but importantly this was 
significantly lateralized to the right hemisphere (Figs 3A and 
5v4). In addition, we dissected a dorsal frontoparietal pathway, 
again more prominent in the right hemisphere, connecting the 
FEF to a region in the inferior parietal lobe (Figs 3B and 5B). 
Finally, we also found evidence for a frontoparietal pathway, 
projecting medially along the cingulum, connecting the SEF to 
the ipsilateral PEF, but this pathway was not demonstrable in all 
individuals (Fig. 36). 

Right hemisphere dominance in the gaze control network 
was found for all individuals (Fig. 4) and confirmed by the group 
average data set for both the SEF to FEF pathway, as well as the 
dorsal frontoparietal pathway connecting the FEF to the inferior 
parietal lobe. Such lateralization in the underlying anatomy 
might explain recent fMRI findings that demonstrate right 
hemisphere dominance in cortical activity in response to visually 
guided saccades (Petit et al. 2009). Consistent with these 
findings, lesions to the right hemisphere following stroke have 
long been known to cause more severe and long-lasting deficits 
in eye movements compared with left hemisphere lesions (De 
Renzi et al. 1982; Perry and Zeki 2000). Indeed, hemispheric 
lateralization has been found for many cognitive functions, most 
notably within the language domain (Powell et al. 2006; Catani 
et al. 2007; Vernooij et al. 2007) and appears to be a uniquely 
defining feature of the human brain — thought to accompany the 
evolutionary advance in cognitive ability (Hopkins and Rilling 
2000; Vallortigara and Rogers 2005). In contrast, anatomical and 
functional aspects of the macaque brain are largely symmetrical, 
and lesions to the right or left hemisphere cause comparable 
deficits (Gaffan and Hornak 1997; Kagan et al. 2010). 

Similar right hemisphere dominance also exists within the 
network for spatial attention in the human brain (Kim et al. 
1999; Mesulam 1999; Bartolomeo et al. 2007; Umarova et al. 
2010) and given the intimate relationship that exists between 
spatial attention and eye movements (Corbetta et al. 1998; 
Nobre et al. 2000), perhaps it is not surprising that we find 
similar lateralization for the oculomotor system. Anatomical 
connections underlying the oculomotor system may well be 
shared with those for orienting attention when an oculomotor 
response is required (Umarova et al. 2010). 

Structural connectivity studies in monkeys have highlighted 
the importance of the long association fibers of the superior 
longitudinal fasciculus (SLF) in spatial processing (Mesulam 
1981; Petrides and Pandya 2006; Schmahmann et al. 2007), 
these fibers connect the frontal lobe to the parietal, temporal 
and occipital lobes, and similar frontoparietal pathways have 
been identified within the human brain (Catani et al. 2002; 
Makris et al. 2005; Thiebaut de Schotten et al. 2005; Doricchi 
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View from below 



Figure 5. Reconstructed tracts for the group average data set. Reconstructed tracts are illustrated in 3D superimposed on coronal, axial, and sagittal slices of the SPM T } 
template image, together with the group fMRI results (RFX analysis, illustrated at P = 0.005), all normalized and coregistered in MNI space. (4) Reconstructed pathways within 
the frontal lobes — tract connecting contralateral SEF shown in red, ipsilateral SEF to FEF connections shown in green. {B) Parietal connections — contralateral PEF connections are 
illustrated in green and the lateral tract connecting ipsilateral FEF and inferior parietal cortex is shown in yellow. Unlike all the other tracts described, this latter tract does not 
terminate in a region of cortex activated by our oculomotor task. Instead this lateral tract terminates in the inferior parietal cortex, near the junction of the angular gyrus and 
supramarginal gyrus (homologous to BAs 39/40): approximate MNI coordinates: L [—38/— 62/46], R [41/-55/46] — an area known to be involved in visually guided saccades. 
This figure was generated using "ExploreDTI" (Leemans et al. 2009). 



et al. 2008; Gharabaghi et al. 2009; Bernal and Altman 2010). 
The right lateralized frontoparietal pathway dissected here, has 
closest homology with the middle segment of the SLF-II, which 
connects the dorsolateral prefrontal cortex to the inferior 
parietal cortex. Damage to this pathway within the right 
hemisphere results in visual neglect and a profound deficit in 
orienting and sustaining spatial attention toward contralateral 
space (Mort, Malhotra, et al. 2003; Bartolomeo et al. 2007; 



Husain and Nachev 2007; Shinoura et al. 2009; Singh-Curry and 
Husain 2009). The right inferior parietal region to which the 
FEF connects is known to be involved in visually guided 
oculomotor tasks (Perry and Zeki 2000; Mort, Perry, et al. 2003; 
Petit et al. 2009; Kagan et al. 2010), when the target for an eye 
movement is visible and needs to be spatially mapped. 
However, this area was not consistently activated by the 
oculomotor task we employed here, which required voluntary 
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eye movements to be performed in the dark, that is, in the 
absence of a visual target. 

Our interpretation relies on the assumption that the output 
of the tractography algorithm is a true reflection of the 
corresponding white matter tracts, and that any reconstructed 
fibers projecting toward a region of cortical activity will 
ultimately synapse in that region and hence are direct evidence 
of an anatomical connection. However, we acknowledge that 
this is a noninvasive technique, and with the currently 
achievable image resolution and the complex architecture in 
juxtacortical white matter, it is possible that fibers that appear 
to project to specific regions of cortex may not actually 
terminate there. Nevertheless, we believe tractography pro- 
vides us with the best possible insights into white matter tracts 
in the living human brain that cannot currently be achieved by 
any other modality. 

Although we have successfully dissected white matter 
pathways connecting key oculomotor regions in the human 
brain, it is of course possible that additional connections exist 
that have not been identified here. For example, a direct 
pathway between the FEF and the superior parietal cortex 
might have been expected, based on anatomical findings in 
macaque (Huerta et al. 1987; Stanton et al. 1995), as well as 
fMRI studies in humans that demonstrate simultaneous activity 
in these regions during eye movements and spatial orienting 
tasks (Corbetta 1998; Natale et al. 2009). It is possible that 
alternative tractography approaches such as spherical decon- 
volution based techniques (Frank 2001; Tournier et al. 2004; 
Dell'acqua et al. 2010; Jeurissen et al. 201 1), may yield different 
results — especially in areas of complex fiber architecture. 
However, such methods have their own limitations and 
shortcomings (for review, see Tournier et al. 2011). For 
instance, bending single fiber bundle populations (i.e., path- 
ways with high curvature at the length scale of the voxel size) 
could be perceived as multiple fiber populations by these 
alternative approaches. As this is an active field of research, it 
currently remains unclear how this will affect the identification 
of specific pathways within complex fiber systems. 

By combining information from both fMRI and diffusion 
MRI-based tractography, we have found strong support for the 
existence of key structural links between cortical areas 
involved in gaze control in humans. This critical information 
not only allows us to explore similarities and differences 
between species, it provides a basis for constraining future 
exploration into the functional role of these connections, and 
whether such connectivity modulates with task requirements. 
Furthermore, within a clinical setting, combined tractography 
and fMRI can be used to assess changes in network 
configuration that might occur following brain damage or as 
a result of subsequent neurological treatment (Grefkes et al. 
2010). Such information could be used to constrain as well as 
monitor the progress of different methods of treatment and 
rehabilitation. 

Previous multimodel imaging studies have defined regions of 
cortex identified using fMRI to constrain tractography algo- 
rithms. However, only one of these studies used an initial 
unconstrained anatomically guided approach, followed by 
fMRI-constrained path reconstruction (Conturo et al. 1999). 
In addition, with one exception (Powell et al. 2006), these 
studies did not seek to establish the existence of hemispheric 
bias, and hence did not control for hemispheric difference in 
the size of the seed ROIs or else chose to examine just one 



hemisphere (Upadhyay et al. 2007; Umarova et al. 2010). Here, 
we performed both anatomically guided and fMRI-constrained 
path reconstructions, using controlled methods to ensure no 
loss of information due to variability in individual anatomy as 
well as preventing the possibility that hemisphere bias in seed 
definition could have been introduced by the operator. 

Conclusions 

We combined diffusion MRI-based tractography with fMRI to 
explore white matter connections within the human oculo- 
motor system. Our data provide compelling support for the 
existence of white matter pathways that have previously only 
been described in detail in nonhuman primates. In addition, we 
highlight an important new finding that appears to be unique to 
humans — that is, a significant right hemisphere bias in the 
white matter pathways connecting oculomotor control cen- 
ters. This lateralization exists for connections within the frontal 
lobes, as well as between the frontal and parietal lobes. 
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